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Many experiments in quantum information aim at creating graph states. Quantifying the pu- 
rity of an experimentally achieved graph state could in principle be accomplished using full-state 
tomography. This method requires a number of measurement settings growing exponentially with 
the number of constituents involved. Thus, full-state tomography becomes experimentally infeasible 
even for a moderate number of qubits. 

In this paper we present a method to estimate the purity of experimentally achieved graph states 
with simple measurements. The observables we consider are the stabilizers of the underlying graph. 
Then, we formulate the problem as: "What is the state with the least purity that is compatible 
with the measurement data?" We solve this problem analytically and compare the obtained bounds 
with results from full-state tomography for simulated data. 



I. INTRODUCTION 

r- 1 ! 

Stabilizer states [l[ (notably including cluster states Q) represent a major class of entangled states. They form 
the resource for a number of applications in quantum information, such as quantum computing [3(, quantum error 
correction [J], and quantum cryptography Q. Because of the usefulness of cluster states, a considerable effort has been 
devoted to their theoretical studies as well as to their experimental implementation. To date, photonic cluster states 
of four [| I and six qubits Q have been experimentally demonstrated, and their realization with trapped ions is actively 
pursued 8| . 

The experimental progress made in the area of cluster states effects the need for sophisticated methods for the 
estimation of the system's properties. A natural problem which arises is the fact that the density matrix grows 
exponentially with the number of qubits involved. So far, it has already been shown that fidelities and entanglement 
in cluster state experiments can be estimated using a number of observables linear in the number of qubits [UEIj. The 
observables of choice are the so-called stabilizers (see e.g. Ref. [TTJ] for an introduction to the stabilizer formalism), 
which are given by 



|>. K j = X J Z Nj) (1) 

O 

where X denotes the usual Pauli x operator acting on qubit j, and Zpf j denotes the Pauli z operator acting on 
neighbours of j defined by the underlying graph. Note that for a graph of n qubits, the set {K\, . . . , K n } generates 
an abelian group, called the stabilizer group. Graph states form the simultaneous eigenvector with eigenvalue +1 to 
these stabilizers. 

In this paper we will show how to estimate the purity tr(p 2 ) in cluster state experiments. We will consider the 
generators of the stabilizer group as the observables. Similar to the estimation of entanglement in Refs. [12| (cfr. Ref. 
[13| for a review of other entanglement estimation methods), the problem is formulated as: "What is the quantum 
state with the lowest purity, and which is compatible with the measurement data?" Mathematically, this question is 
the following quadratic optimization problem: 

Pram = min p [tr(p 2 ) : tr(Kip) = a*, p > 0] . (2) 

The paper is structured as follows: in Sec. Ill Al we will utilize the symmetries allowed by the observables to restrict the 
optimization to stabilizer diagonal states. Then, we provide an exact analytical solution to the above optimization 
problem in Sec. Ill Bt The following section proceeds with a discussion of the quality of the obtained solution by 
comparing it with results from full-state tomography from simulated data of a noisy system. Sec. IVII discusses the 
estimation of the von-Neumann entropy from stabilizer measurements, and provides an analytical lower bound on the 
entropy. 
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II. ESTIMATING PURITIES 



A. Symmetries of the stabilizers 



As mentioned in the introduction, the natural observables in cluster state experiments are given by the stabilizers. 
Let us assume the goal of an experiment was the creation of a pure cluster state, with stabilizers Kj, j S {1, . . . , n}. 
We denote the measurement outcomes by a,j = tr(pKj). Note that these expectation values are invariant under 
rotations of the stabilizer group. By twirling over the stabilizer, one may therefore restrict to states of the form 

P=^ E '•<:...< /V... A',', • (3) 
ii,...,i„=0 

Here, the following twirling protocol was utilized: 

p'^P=^ E K i ?...K%p'K?...K i j. (4) 

ii,...,i n =0 

The coefficients in the stabilizer decomposition of p are subject to the normalization constraint C0...0 = 1 and measure- 
ment outcomes a 3 = tr(pKj) — c^, where j denotes a bit string of zeros, but 1 in position j. Furthermore, any valid 
density matrix must be positive-semidefinite. The stabilizer decomposition of p provides a convenient way to calculate 
the eigenvalues of p. Since the stabilizers mutually commute, and the spectrum of each term in the decomposition is 
simply given by {+1,-1}, one finds the following expression for the eigenvalues: 



E M) E < 



H Jl, 



(5) 



B. Lowest purity compatible with measurement data 



In this section we determine the lowest purity compatible with measurement data in cluster state experiments. So 
we calculate the exact solution to the optimization problem @. As seen in the previous section, it is legitimate to 
restrict to stabilizer diagonal states of the form Q . Then, the problem boils down to solving the following positive 
quadratic program: 

I 1 

minimize P(c) = — E c n-t n ( 6 ) 

ii,...,i„=0 



s.t. cj = dj (7) 

-\, ' (8) 

Clearly, this problem can be solved numerically by convex and quadratic optimization solvers such as SDTP, Sedumi, 
or Yalmip. However, the number of coefficients Cjj...j n grows exponentially with the number of qubits. Therefore, a 
numerical approach becomes practically intractable even for moderate qubit numbers. We will therefore derive an 
analytical solution. 

Without loss of generality, we restrict the measurement values to be positive. Furthermore, we will only consider 
the case, in which the measurement outcomes of neighboring qubits fulfil + ft& 2 > 1. In other words, we restrict 
to the case of a non-zero fidelity with the desired stabilizer state. 

Now one can easily check that the choice 

<Hl...i n = E ikClk ~ E ik + 1 ( 9 ) 

k k 

fulfils all the constraints in the above quadratic program. Positivity of the eigenvalues follows from the fact that p is 
diagonal in the stabilizer basis. Thus, these coefficients represent a solution, which we denote by c* in the following, 
to the primal problem. 
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Since the primal problem is convex, we may simply check the Karush-Kuhn- Tucker (KKT) conditions for optimality 
of the solution |l4|. To be explicit, the KKT conditions require 

A, ....,•:<-:. <>. (10) 
g(c 3 = c|) := c- - a, = 0, (11) 

/', > 0, (12) 

Vjl...j n ^jl...jn( C *) = (!3) 

1 n 

Vh-jnP(c*)~ E V / ./ A < • ' E ./ •■>"■ '• ( 14 ) 

ii,...,i 7 i=0 m— 1 

The first two conditions state that c* is primal feasible, which is already proved. The last condition guarantees the 
optimality of c*. Instead of deriving the dual quadratic program, we will now prove that Lagrange multipliers 
and Vm can always be found to the solution c*, such that the KKT conditions are fulfilled, thus showing that c* is 
indeed optimal. 

A valid solution for the Lagrange multipliers is provided in the following way: in order to satisfy (|13[) choose 

Hi-jn = °> if X ji-jA c *) > 0- 

Condition (p~4|) can be rewritten as 

2 1. 

— c-— Au + v = 0, (15) 
2n 2 n 



where it follows from Eq. [5] that A is the normalized Hadamard matrix with elements A ix ^ n j x m j n — (— 1)^*<= **■'* , 
c = {cj 1 ...j Tl ), A = (Ajj . j^), and v is the vector containing Vf n for m = 1,. .. ,n, while all other entries Vj 1 ...j Tl = 0. 
Here the matrix A is chosen such that Ac — 2"A. Because of the orthogonality relation of the Hadamard matrix, and 
of A = A T , it follows with A T A = 2 n I: 

2 

/i = A T (— c + v) =2\ + Av. (16) 

If Aj > (j = (ji, . . . ,j n ))i then it must hold that /ij = 0. Especially, in any case /x^ = 0, since the largest eigenvalue 
of p is Aj. Therefore, we have: 

T£(-l)r* UA = -2\ I . (17) 

Here, / denotes the set of bit-strings with at most one bit unequal zero, as in the vector v only the first n entries are 
unequal zero. More formally, I = {(0, . . . , 0), (1, 0, . . . , 0), . . . , (0, . . . , 0, 1)}. This system of equations has the solution 

v t = \ t - A 5 , j> 0, (18) 

^ = -2A --5>f. (19) 

It remains to show that /ij > for Aj = 0: 

^ = E(- 1 ) 7 '™^ = v ° + E (~ 1 ) f ™^ (2°) 

™ rh=i0 

= -2A ff -^^+^(-l)^ (21) 
= -2A Q -+^((-ir»-l)^ (22) 

m^O 

= -2A 3 + (Mr™ - - Ag) (23) 



m#0 

-2A -+1)^((4)^ - 1) - % (24) 

- 2A o + ^E (-^i J-™) E(- 2(5 M-) C ?- ( 25 ) 



m^0 
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Since it is clear that A<j > 0, we can restrict our attention to eigenvalues Aj = and dual variables pij with at least 
one index j a = 1. Therefore, we have 



E(- 25 u-™)E(- 25 i,?-™)^ c ^ 4 E ^,., Q _ lllQ+ ,.. 4 „ (26) 

m^O i ~ i^=0,/9e{l,...,n}— a 

The RHS of the last inequality can be rewritten as 

4 1 1 1 

2n ^ ] Cii...i Q _ili Q + i...i„ — 4(Ag — ^ ] c ii...i Q1 _i0i a+ i...i„)' (27) 

t^5=0,/9e{l, ...,n}-a t^=0,/3e{l, ...,n}-a 

Employing Eq. the last term may be evaluated as 



1 1 1 1 " 

^ 51 ci 1 ...i a _ lWo+1 .„ iB = 51 ^r(E ifcafc -E* fc + 1 ) ( 28 ) 

i =O,0e{l,—,n}— a i/3=0,l3e{l,...,n}-a k=l k 

1 " A- 

= I (E a fc-" + 2 )<f • ( 29 ) 



4* 



The last inequality follows from the fact that Ag = §(X}fc=i flfc — + 2). Applying this result to Eq. (|2l))) immediately 
implies /ij > 0, if Aj = 0. Therefore, all KKT conditions are fulfilled, proving that the solution c* is indeed optimal. 



III. EXAMPLE: TWO QUBIT PURITY 

For the purpose of illustration, let us consider the simple example of two qubits, supposedly prepared as a cluster 
state. Given only the outcomes of the stabilizer measurements a\ — tr(pX <g) Z) and a 2 — tr(pZ <S> X), we seek the 
lowest purity compatible with these two measurements. Using the symmetries of the stabilizers, we may restrict the 
problem to density matrices of the form 

P = ^(cool + c w Ki + c Q iK 2 + cuKiK 2 ). (30) 

Because of tr(p) — 1, c 00 = 1. Furthermore, the measurements determine Cio = &i and Coi = a 2 . In order to find the 
minimal purity, the objective is now to minimize 

Prain = ^( C 00 + C 10 + C 01 + C ll) (31) 

subject to A JU2 (c) = E! 1 ,,: 2 =oM) 1iJi+12J2 cm 2 > 0. 

As solution to the primal problem, one chooses en = \{a,\ + a 2 — 1). Then the eigenvalues of p are simply given by 

Aoo = ^(«i + a 2 ) (32) 
Aio = ^(l-a 1 ) (33) 

Aoi = i(l-a 2 ) (34) 

An = (35) 

thus fulfilling the positivity constraint. To see the optimality of the solution, one can check the remaining KKT 
conditions. 

- - Mi - M2 - M3 - M4 + v x =0, (36) 

— - Hi + i-p 2 - M3 + M4 + v 2 =0, (37) 

— - hi - [i 2 + p 3 + fi4 + vj, =0, (38) 

i(ai + a 2 - 1) -Mi + M2 +M3 - m =0 (39) 
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TABLE I: Comparison between exact and estimated purities for noisy cluster states 



No. qubits 


exact purity 


estimated purity 


relative deviation 


2 


0.8269 


0.8233 


0.0044 


3 


0.7520 


0.7417 


0.0137 


4 


0.6838 


0.6646 


0.0281 



The Lagrange multipliers v are determined by 

vva = Aio - A 00 = ~(1 - 2ai - a 2 ), (40) 

Vol = Aqi - A 00 = ~(1 - 2a 2 - ai). (41) 

Voo = -2A 00 - uio - vq X = + a 2) - 1. (42) 

If the measurement outcomes fulfil &i,a2 > 0, which should be the case in a cluster state experiment, then 
Aoo,Aio,Aoi > 0, and An = 0. Therefore, pi = p 2 = /Z3 = 0, to fulfill the KKT condition p3[) . It is now 
straightforward to check that the choice of these Lagrange multipliers satisfy conditions (J2SJ) through ([55]). thus 
proving optimality of the solution. 

IV. QUALITY OF THE BOUNDS 

In this section we will test the usefulness of the obtained bound on the purity. For this purpose, let us consider a 
perfect cluster state, which is subject to dephasing for a certain time, which is here assumed to be 10 ms. The system 
is then described by a master equation of the form 

p = \YS- Zi P z *-p)> ( 43 ) 

i 

where 7 is the dephasing rate, which is assumed to be (10 ms) . 

V. ESTIMATION WITH ERROR BARS 

Needless to say, due to a finite number of measurements and experimental imperfections, the measurement outcomes 
ctj we considered in previous sections posses an error Aa^ > 0, which clearly affects the result of a purity estimation. 
In this section we briefly address the question how to take into account such errors. 

Let us assume we measured the stabilizers K% with measurement outcomes a% ± Aai = tr(pKi). As seen in Sec. 
IIIB1 the optimal solution to the purity minimization problem @ is given by a stabilizer diagonal representation of 
the density matrix with coefficients c% x ...i n — ^2 k ikO-k — + From this, it can be easily seen that the choice 

c it...i = J2k^k(ak ± Aafc) — X)fe*fe + 1 leads to upper (lower) error estimates of the minimal purity with objective 

values P^ = ^El 1 ,...,i n =o(4 ±) . i J 2 " 

VI. ENTROPY ESTIMATION 

Alongside the purity of a quantum state tr(p 2 ), the entropy is another quantifier for the degree of disorder of a 
system. Given only the measurement outcomes of the stabilizer, one can employ similar techniques as in the case of 
the purity to estimate the entropy. 

The von Neumann entropy is given by S = —tr(plogp). Again, let us consider the generators of the stabilizer group 
as our observables with measurement outcomes a\ = tr(pKi), . . . , a n = tr{pK n ). Then we may restrict to states 
diagonal in the stabilizer basis p = ^ Ci 1 ...i n K\ 1 ...K^ with eigenvalues Xj 1 ...j n (p) = ^2n x j =o( — ■"•) h **^* c ii...* n ' 
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TABLE II: Comparison between exact and estimated entropies from the least purity state for noisy cluster states 



No. qubits 


exact entropy 


estimated entropy 


relative deviation 


2 


0.3827 


0.3803 


0.0063 


3 


0.5740 


0.5667 


0.0127 


4 


0.7653 


0.7505 


0.0193 



Now, we need to maximize the entropy over all states predicting the measurement outcomes Ofc, k € 1, ...,n. 

S m ax =max{S(p) : tr(pK t ) = a { ,p> 0} (44) 
l 

=min{ 2J A, lof/X, : A^...^ > 0, cj = a 3 } (45) 

i±,...,i n = 

Clearly, the solution c^...^ = ^2 ik<ik — ^2 ik + 1 satisfies the constraints, thus providing a lower bound on the entropy. 
Considering initially perfect cluster states subject to dephasing according to Eq. ([43|) . one obtains the results given in 
the table ITT] Remarkably, the estimate differs only a few per cent from the exact value of the von Neumann entropy. 

Even though the least purity state provides a very good lower bound on the maximal entropy, an exact derivation of 
the entropy is possible in the following way: given the probabilities pj^ = 1 ^ afc , then it follows from the subadditivity 
of the entropy S({pk}) < S(pk), that the maximal entropy is given by 

l 

Smax = ~ A...... lof/X, . (46) 

ii,...,i n =0 

with Ajj.,.1,, = rifci a i an d the indices au being a binary index. 



VII. CONCLUSION 



In this paper we have developed techniques to estimate purities and entropies in cluster state experiments with few 
measurements. We assume only that the generators of the stabilizer group are measured in an experiment. Therefore, 
the number of local measurements is linear in the qubit number, compared to an exponential growth in the number 
of measurement settings for full-state tomography. Given only these measurement outcomes, we have developed the 
optimal analytical solution to the minimal purity compatible with the measurement data. A comparison of this bound 
on the purity with the exact purity calculated from simulated data of noisy cluster states shows that, despite the little 
number of measurements, the bound on the purity differs only by a few per cent from the exact value. In a similar 
way, we obtained a lower bound on the von Neumann entropy considering only measurements of the stabilizers of a 
cluster state. This bound deviates only by a few percent from the exact value of the entropy, too. In addition, we 
employed the subadditivity of the entropy to obtain an analytical result for the exact entropy of a stabilizer state. 
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